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ABSTRACT 

Aims. Most of the known debris discs exhibit cool dust in regions analogous to the Edgeworth-Kuiper Belt. However, a rare subset 
show hot excess from within a few AU, which moreover is often inferred to be transient from models for planetesimal belt evolution. 
In this paper we examine 2 such sources to place limits on their location to help distinguish between different interpretations for their 
origin. 

Methods. We use MIDI on the VLTI to observe the debris discs around q Corvi and HD69830 using baseline lengths from 44— 130m. 
New VISIR observations of HD69830 at 18.7/im are also presented. These observations are compared with disc models to place limits 
on disc size. 

Results. The visibility functions measured with MIDI for both sources show significant variation with wavelength across 8- 13pm in 
a manner consistent with the disc flux being well resolved, notably with a dip at 10-1 1.5/im due to the silicate emission feature. The 
average ratio of visibilities measured between 10-1 1.5/jm and 8-9/jm is 0.934+0.015 for HD69830 and 0.880±0.013 for q Corvi over 
the four baselines for each source, a departure of 4 and 9<r from that expected if the discs were unresolved. HD69830 is unresolved by 
VISIR at 18.7/jm. The combined limits from MIDI and 8m imaging constrain the warm dust to lie within 0.05-2.4AU for HD69830 
and 0. 1 6-2.98AU for q Corvi. 

Conclusions. These results represent the first resolution in the mid-infrared of dust around main sequence stars. The constraints 
placed on the location of the dust are consistent with radii predicted by SED modelling (1.0AU for HD69830 and 1.7AU for q Corvi). 
Tentative evidence for a common position angle for the dust at 1.7AU with that at 150AU around q Corvi, which might be expected 
if the hot dust is fed from the outer disc, demonstrates the potential of this technique for constraining the origin of the dust and more 
generally for the study of dust in the terrestrial regions of main sequence stars. 

Key words, circumstellar discs - main-sequence stars 



1. Introduction 

Debris discs are believed to be remnants of planet formation. 
As well as providing a unique window on how a system may 
have formed and evolved, studies of debris discs can reveal much 
about the current status of a system: the location of the dust can 
indicate possibl e configuratio ns of any planetary system around 
the star (see e.g. Wyatt 2008); and the amount of debris material 
allows us to understand what physical conditions ma y be like 
for any habitable planets in the system (Greaves 2006). Debris 
disc emission is seen to peak typically longwards of 60/im im- 
plying that the dust is cool <80K and thus predicted to lie at 
regions many 10s AU from the star. In the majority of cases 
where this emission has been resolved around Sun-like stars 
it has been confirmed to l ie at >40AU (see e.g.. iHolland et alj 
119981: iGreaves et alj|2005 l : iKalas et alj|2007l) . Exoplanet studies 
on the other hand typically probe close-in regions for giant plan- 
ets, as most planets detected to date have been found through ra- 
dial velocity technique s with their inherent detection biases (see 
e.g. iButlereTaT] 120061) . Debris disc studies are therefore typi- 
cally complementary to exoplanet studies. Furthermore, com- 
parison with the Solar System's debris discs, the asteroid and 
Edgeworth-Kuiper belts, can help us further understand the like- 



lihood of finding Solar System analogues in our solar neighbour- 
hood. 

Very few stars exhibit hot dust within ~10 AU, i.e. in the 
regio n where we expect planets may have formed (llda & Linl 
2004). Four surveys have searched for hot dust around Sun-like 
stars by looking for stars wit h a 25 tan flu x in ex cess of photo- 
spher ic levels usin g IRAS dGaidoslll999|). ISO (|Laureijs et all 
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2002) and Spitzer Mines et alj 120061: iBrvden et al.ll2006l) . All 
concluded that only 2 + 2% of Sun-like stars have hot dust with 
infrared luminosities / = L/«/L* > 10~ 4 . These hot dust sources 
provide a unique opport unity to prob e the terrestrial plan et re- 
gion of their systems (e.g. lWvattll2002l) . IWvatt et al.ld2007l) iden- 
tified 7 main-sequence Sun-like (F, G and K-type) sources with 
hot dust emission at ~1AU. For these sources the levels of ex- 
cess emission observed were compared with the predictions of 
analytical models for the collisional evolution of debris discs. 
This model predicts that there is a maximum level of emission 
that a disc of a given age and radius can have; discs that are 
initially more massive, which one may assume might end up 
brighter, process their mass through the cascade more quickly 
(i.e. more massive di s cs hav e shorter collisional timescales). 
Although lLohne et ail J2008) found that the analytical model 
may underestimate the levels of dust emission produced by de- 
bris discs at ag es of ^ 1 Gyr by a fac tor of 10, they concurred with 
the findings of Wy att et all ((2007) that for 4 of these 7 sources 
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the level of dust emission is too high to be explained as the 
product of a steady-state collisional cascade in a planetesimal 
belt coincident with the dust at ~1AU, and must therefore be the 
product of a transient event. Recent near-infrared interferometric 
observations have also revealed 4 sy stems with hot d ust emis- 
sion at -0.1 AU fro m the star (Vega, lAbsil et al.ll2006t t C eti. 
IdiFolco et af] 120071: /? Leo and £ Lep. lAkeson etal.ll2009l al- 
though the hot emission around ( Lep must be confirmed). Dust 
lifetimes in these regions are very short and these dust popula- 
tions are also likely to have been produced in a recent transient 
event. 

A critical diagnostic to understand the origin of this transient 
emission is the location of the disc, and in particular whether the 
emission is asymmetric. In this paper we examine two sources 
inferred to have transient emission at ~1AU, HD69830, a K0V- 
type star with an age of 2Gyr (Beichma n et al.l 120061) and n 
Corvi , an F2 V-type star with an age of 1 .3Gyr dMallik et al.l 
120031) . The predicted dust locations around both sources is 80- 
90 mas (from SEP modelling of the IR S spectra of the targets, 
Beichman et al. 2005t lChen et al.f 2006). although there remains 
considerable uncertainty in the true radial location of the emis- 
sion. Thus far 8m imagi ng of 77 Corvi sug gests the dust lies 
within 3AU of the star dSmith et al.l 120081) . Here we present 
VISIR observations of HD69830 and MIDI observations of both 
targets. 

This paper is structured as follows. In section 2 we describe 
our interferometric observations and data reduction procedures. 
In section 3 we describe the results of the MIDI observations. In 
section 4 we present new VISIR imaging of HD69830 showing 
that the excess emission at 18.72/vm is not-extended, and discuss 
the implications for the maximum extent of the disc. In section 5 
we discuss the results and the implications for the radial location 
of the debris discs. We conclude in section 6. 



2. MIDI observations and data reduction 

For readers unfamiliar with the basics of optical/infrared inter- 
ferometry, an outline of the basics and more detailed information 
on the VLTI in particular can be found i n Observation and Dat a 
Reduction with the VLT Interferometer, (M albet & Perrinll2007l) . 



2.1. Observations 

Interferometric data were obtained on HD69830, rj Corvi and 
suitable standard stars in visitor mode at the VLTI under pro- 
posal 078.D-0808 on 5th-9th March 2007. The expected levels 
of flux from our target sources in the N band (see Table Q] and 
Figure Q) required the use of the 8m UTs and the HIGH-SENS 
observing mode (for details see below). As far as was possible 
we tried to observe on both short and long baselines at near per- 
pendicular position angles so as to constrain the size and geome- 
try of the sources. In addition, observations of the science targets 
were bracketed with observations of standard stars to assess any 
variability of the interferometric transfer function. The science 
and standard star observations in this program are summarised 
in Table [2] 

In all cases data were secured following the sequence of 
steps outlined below. A more detai led discu s sion o f this, together 
with its rationale can be found in iTristraml (I2007I) . to which the 
interested reader is referred for greater detail. 

- Initial acquisition, for pointing and beam overlap optimiza- 
tion, was performed in imaging mode using the N8.7 filter 
(central wavelength 8.64 fim, width 1.54 /im). We used a 



Table 2. Observations summary 



Date and 


Baseline 


Pos. Angle 


Seeing 


Source 


Target 


Config. 


m 


EoN 




HD 


type 


05/03/2007 


98.9 


40.5 


0.60 


61935 


Cal" 


T IT 1 T ITO 

UT1-UT3 


96.4 


40.9 


0.55 


69830 


Sci 


U I 


1 no 1 




n qc 


yjZ/Z 


cal 




1 n 1 


39.5 


1 in 
1.5V 


rj Corvi 


Sci 




102.4 


38.0 


1.50 


116870 


Cal 


07/03/2007 


59.5 


113.4 


1.15 


61935 


Cal 


T TTIO T IT A 

UT3-UT4 


60.4 


113.7 


1.30 


69830 


Sci 


A 


59.8 


114.3 


1.25 


73840 


Cal 6 




n 


1U5.4 


n 

U. /5 


II08/U 


Cal 




62.4 


109.0 


0.80 


r] Corvi 


Sci 




62.1 


112.5 


0.65 


107218 


Cal 


08/03/2007 


130.2 


63.1 


0.90 


61935 


Cal 


UT1-UT4 


130.2 


62.8 


0.75 


69830 


Sci 


C 


130.1 


63.3 


0.75 


73840 


Cal 


TTT1 TTT^l 
U 1 1-U Lj 


1 nn 1 


jV.J 


1 1 c 


1 loo /U 


t-.ai 


D2 


102.3 


38.3 


1.60 


t] Corvi 


Sci 




102.2 


36.5 


0.70 


116870 


Cal 


09/03/2007 


42.7 


32.5 


1.50 


61935 


Cal 


UT2-UT3 


44.4 


43.8 


1.55 


69830 


Sci 


B 


44.2 


34.3 


1.20 


73840 


Cal 




42.2 


26.8 


0.65 


116870 


Cal 




46.1 


39.2 


0.75 


t] Corvi 


Sci 




46.6 


42.4 


1.15 


107218 


Cal 



Baselines and position angle (projected on the sky) are given as the 
mean for the observation. 

" The first standard star observation on run Dl (of HD61935) showed 
great variation in visibility with wavelength (see section 3.4), and so 
this was rejected as a suitable observation for calibration. 
* During run A the observation of standard star HD73840 could not be 
used, as the telescope beams showed poor alignment on the readout ar- 
ray. Thus the observation of HD69830 was calibrated with HD61935 
only, since the next usable observation of a standard (HD 116870) was 
taken much later in the night after sky conditions had changed signifi- 
cantly. 



standard chop throw of 15"for acquisition and observations 
of the total source intensity. 

Subsequently, fringe searching was performed using unfil- 
tered but dispersed light. A slit width of 0V52 width was used 
together with the MIDI prism dispersing element giving a 
spectral resolution, A/AA, of 30 at A — 10.6pm. 
After fringe centering, longer sequences of fringe data were 
secured. For the first few frames of each sequence, the de- 
lay line positions were altered so as to be significantly offset 
from the zero optical path difference (OPD) position to al- 
low a determination of thermal and incoherent backgrounds. 
Thereafter, the delay lines were returned to their nominal lo- 
cations and typically a few hundred scans of fringe data were 
secured. These data allowed an ongoing determination of the 
OPD corrections needed to compensate for uncorrected side- 
real and atmospheric motions. It this way it was ensured that 
that the fringe packets always remained well centred. 
Finally, at the end of each interferometric data sequence, the 
same instrument set-up was used to determine the through- 
put of the individual UT optical trains. For these "intensity 
calibration" measurements chopping of the secondary mir- 
ror was used to suppress any thermal background. The chop 
frequency and integration time before readout were varied 
between the runs to try to maximise the signal-to-noise on 
these measurement on a case by case basis. 
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Science targets 


Source 


Spec type 


Age 


RA 


Dec 


F, 


at lOyum 


F disc at lOyum 


Predicted disc size 


HD 




Gyr 








mJy 


mJy 


mas 


69830 


KOV 


2" 


08 18 23.95 


-21 37 55.8 




872 


102 


80 


r] Corvi 


F2V 


1.3* 


12 32 04.23 


-16 11 45.6 




1736 


371 


90 


Calibrators 


Source 


Spectral type 


RA 


Dec 




F, at 10 yum 


Angular size 


HD 












mJy 


mas 


61935 


G9III 




07 41 14.83 


-09 33 04.10 




9490 


2.24±0.01 


73840 


K3III 




08 40 01.47 


-12 28 31.30 




12312 


2.40±0.01 


95272 


Kim 




10 59 46.46 


-18 17 55.56 




9510 


2.24±0.01 


116870 


K5III 




13 26 43.17 


-12 42 27.60 




10416 


2.58±0.01 


107218 


M4III 


12 19 42.59 


-19 11 55.97 




7913 


2.16±0.02 



For science targets the stellar flux was determined by Kurucz model profiles scaled to the 2MASS K band photometry. Total fluxes taken from 
IRS photometry of the targets were then used to determine the disc flux, with the predicted size being based on SED fitting with a blackbody. 
For calibration t argets the angular size is as given by the CalVin tool available at http://www.eso.org/instruments/midi/tools Key: " Age from 
Beic hman et al.l J2006h : b Age from lMallik et alJd2003h: note that the X-ray luminosity of this star is close to the mean value for the Hyades cluster 
which may suggest a younger age of 600-800Myr dStern et alj 19951) . 




Wavelength, fj,m Wavelength, fj,m 

Fig. 1. The emission spectra of o ur science targets in the MIDI wavelength range. Left: The Spitzer IRS spectrum of HD69830 
presented in Beich man et alj d2005l) shows that the excess emission has a large silicate feature in the N band. Right: The Spitzer IRS 
spectrum of 77 Corvi also shows a large silicate feature at ~ 1 1/j.m. This spectrum was originally presented in lChen et al.l (12006). For 
both sources the emission expected from the stellar photosphere is taken from a scaled Kurucz model profile. 



In summary, each full dataset for a single observation of a 
target comprised a fringe file (containing the interferometric data 
from the combined beams) and 2 "intensity calibration" obser- 
vations (one for each beam). Similar data were also secured on 
a standard star (i.e. of known flux and small diameter) so that 
flux calibrated visibilities could extracted. The standards were 
selected from the spectro-photometric catalogue available on- 
line at http://www.eso.org/instruments/visir/tools and are listed 
in Table □ 



gave a poor fit to the peak of the detected emission and were 
much broader than the emission in the direction perpendicular 
to the spectral dispersion. We therefore used masks determined 
by a fit to the total intensity frames as provided by the MIA re- 
duction package instead (see Figure We then used the EWS 
group delay analysis to align the fringes before vector averaging 
over time to derive the correlated flux (for details see the EWS 
user manual). The correlated flux 7 C01T was then compared to the 
total source flux, 7 tot , to give the source visibility V = 7 C orr Aot- 



2.2. Reduction with the MIA+EWS package 

Data reduction was performed using the EWS soft- 
ware, available as part of the MIA+EWS package (see, 
http://www.strw.leidenuniv.nl/~nevec/MIDI/index.html and the 
manual and details therein). 

The first step in the data reduction involved compressing the 
fringe and total intensity frames in the direction perpendicular 
to the spectral dispersion to obtain one-dimensional fringe and 
intensity spectra. Prior to compression, each frame was multi- 
plied by a mask to reduce the impact of noise. We found that the 
default EWS masks provided for HIGH-SENS mode observing 



2.3. Total Source Intensity 

Total intensity data were obtained with MIDI after the fringe 
observations had been secured as described in Sec. 2.1. The re- 
sults showed a great degree of variation between different ob- 
servations of the same object, primarily because of background 
which could not be perfectly subtracted. We found the back- 
ground did not vary linearly in the direction perpendicular to the 
spectrally-dispersed direction, and moreover varied with time 
during the observations, suggesting that sky (and/or telescope) 
variations were significant. As the total intensity measurements 
of the same science target could vary by >50%, we chose to use 
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Fig. 2. Total intensity from MIDI observations. Left: Coadded intensity frames from an observation of a standard star in a single 
optical train (in this case HD61935 run Dl beam "A" = telescope UT3). The spectral dispersion runs along the horizontal direction 
(increasing wavelength from right to left). Right: A cut through one of the spectral channels (at 8.5/im). This shows the profile of 
the total intensity along the pixel direction. Overplotted are the EWS default mask (fixed based on observation set-up) and the MIA 
mask (determined from the intensity image itself). It is evident that the MIA mask fits not only the width of the true detected profile 
but also the location of the peak on the array much more accurately than the default mask. 



available IRS spectra of our science targets as the total intensity 
measurements in ste ad. The IRS spectra (o riginally presented in 
IChen et al . 2006 and Beich man et al.l2005l) were extracted along 
a slit of width 3'.'7 or 4'.'7 depending on the spectral resolution. 
Our MIDI data were secured with a slit width of 0752 and so it 
was necessary to check whether the Spitzer aperture might have 
included additional photometric signals not seen by MIDI. 

For both sources we found there were no background or 
companion objects in t he Spitzer aperture that might be excluded 
by the MIDI aperture dSmith et al.ll2008l) . Th e 3 Neptune mass 
planets discovered within 1AU of HD69830 (lLovis et al.i r2006) 
are not expected to b e brighter in the mid-infrared than 10~ 3 F* 
dBurrows et al J 120041) . and so would only contribute a discrep- 
ancy of at most ~lmJy at 10/vm to our tot al intensity mea sure- 
ment. A full analysis of the IRS spectra dLisse et al]|2007h re- 
vealed evidence for water ices which, due to their lower temper- 
atures, may reside at larger radial offsets from the star than the 
bulk of the emission around HD69830. However, as pointed out 
by the authors the local thermal equilibrium temperature for the 
region of their best fitting model — a disc at ~ 1AU — is only 
245K, cool enough that should the water ice be isolated from the 
hotter dust particles it could be stabilised by evaporative sub- 
limation. In addition a wide Kuiper belt-like l ocation is ruled 
out by the 70yum limits on the excess (l±3mJy, Beichman et al. 
2005). The contribution of the water ice in the 8-13/mi region 
of the spectrum is at the level of 1-3% of the excess emission 
between 8-1 l//m, rising to 27% of the excess emission at 13 /im. 
Thus we believe the effect of the water ice contribution falling 
outside the MIDI beam can only have a significant effect on the 
measured flux longwards of 11.5^m, if at all. The temperature 
of the dust around rj Corvi is predicted by SED fitting to be at 
around 320K, suggesting a small radial offset from t he star . The 
alternative fit to the SED suggested by IChen et alj (|2006) sug- 
gests there may be dust at two temperatures of 360K and 120K, 
but the 8-13/im range is dominated in this model by amorphous 
olivines at the hotter temperature. Thus in either the single tem- 
perature fit or the fit proposed by IChen et al.l (|2006) the excess 
emission in the MIDI wavelength range is expected to be domi- 
nated by grains close to the star. Therefore for both science tar- 
gets it is unlikely that significant emission appears in the IRS 
spectra that would be excluded from our MIDI observations. 



2.4. Visibility calibration 

The total source intensities of the standards, which were chosen 
for their lack of infrared spectral features as well as their small 
angular size, were modelled by Rayleigh-Jeans functions scaled 
to the 10.5^/m flux listed in the spectro-photometric catalogue 
(see Table [T). The visibilities of the standard stars were calcu- 
lated assuming these could be modelled as uniform discs with 
diameters as listed in Table Q] 

The calibration of the correlated flux and determination of 
visibilities for each observation were determined as 

^con^tar A:orr,tar/^coiT,cal ^ ^tot.cal ^ ^cab (-0 

where 7 C on-,tar and / corr , ca i are the correlated fluxes in ADU (output 
of data reduction procedure) of the target and standard star used 
as a calibrator respectively, F conMr is the calibrated correlated 
flux of the target, F to t, C ai is the known flux of the standard star 
(from Rayleigh-Jeans slope) and V ca i is the visibility of the stan- 
dard star which is dependent on the baseline length and wave- 
length. The visibility of the target is then 

V — ^corr,tar/^*tot,tar> (2) 

(Ftot tar again from Rayleigh-Jeans slope). For the case of the sci- 
ence targets F tol t . dr was given by the IRS spectrum. We used the 
average of the two standards bracketing each science observation 
to determine the calibrated visibilities, and difference between 
the two to determine our absolute calibration uncertainties. 

We also performed the above calibrations with the standard 
star observations as "dummy targets" to: 1) determine the level 
of accuracy to which we could trust the absolute visibility mea- 
surements; and 2) examine the visibility functions for evidence 
of changing visibility with wavelength. Examples of the transfer 
function of standard - standard calibration (where the standards 
were observed close in time and space to each other) are given in 
Figure[3] In the left panel we show the transfer function; the vis- 
ibility determined for our "dummy targets" as a function of the 
visibility predicted for this target ( V ca i in equation[T]i. The abso- 
lute levels of this transfer function were found to have a mean 
and standard deviation of 0.99 ±0.10 (Figure [3] left) after cal- 
culating the weighted mean of each visibility over the 8-13jt/m 
range (excluding the 9.2-10//m region which is dominated by 
ozone emission). The weights were determined by splitting the 
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8 9 10 11 12 13 8 9 10 11 12 13 

Wavelength, fim Wavelength, fim 



Fig. 3. Plots of visibilities determined for standard star targets using other standard stars. Pairs of standard stars observed closely 
in time were considered as these are the pairings used for science observation calibration. Note that the shaded region 9.2-lOpm 
is prone to strong uncertainty due to ozone absorption, and so should be ignored. For each pair the 'target' and 'calibrator' were 
chosen randomly. Left: The calibrator-calibrator determined visibilities divided by the predicted visibility for the target star. Perfectly 
calibrated data would give a straight line at unity with no variation with wavelength. We used the weighted mean (over wavelength) 
for each calibrator-calibrator dataset to quantify the quality of the calibration and found a range between 0.86 and 1.12. Right: 
Checking the absolute variation with wavelength for calibrator-calibrator visibilities. All visibilities were scaled to 1 in the range 8- 
9fim. The ratio of visibility between 10-1 1.5pm and 8-9/mi was calculated to determine the variation of visibility with wavelength. 
The maximum deviation from unity for this value was 0.047, with some calibrations providing a much flatter visibility. 



fringe data into 5 equal length sets and calculating the variance 
between these subsets of data at every spectral channel. This ab- 
solute variation is the same as the 10% error typically expected 
in ca libration due to changing conditions between observations 
(e.g. lChesneaul2007h . The visibility functions were fairly flat, as 
shown in the right-hand panel of Figure [3] The weighted mean 
of each visibility function between 8-9/im was used to scale the 
visibilities to 1 in the 8-9/mi range for this plot. The weighted 
mean for each function between 10-1 1.5pm and 8-9pm were 
then compared. This ratio was found to be have a mean and stan- 
dard deviation of 0.991 ± 0.037. Thus the differential visibility 
(i.e. the visibility as referenced to that at some fixed wavelength) 
is typically known with ~ 3 times better accuracy than the abso- 
lute visibilities. 



will vary solely on the ratio of stellar to total flux whereas the 
second will depend both on the ratio of disc to total flux and 
on whether the disc is resolved by the interferometric baseline. 
More precisely, the visibility function will tend to the value of 

Vto, I/ _ ' V„ = ^- (4) 

^disc^O r tot 

when the disc is fully resolved. As can be seen in Figure Q] the 
value of y, es = y*- is expected to change across the spectral 
range of MIDI, particularly in respect of the silicate features seen 
in the excess emission around both stars. As a result, any varia- 
tion in observed visibility as a function of wavelength will need 
to be corrected for this behaviour before it can be interpreted in 
terms of any resolved source structure. 



3. MIDI results 

We now present the results for each of the science targets in turn. 
To help interpret our results, we consider how the visibility func- 
tion of a source is expected to behave if the source is made up 
of multiple components. HD69830 and rj Corvi have significant 
stellar emission in the N band (see Table [TJ. Thus there will be 
two components to the visibility function, the component from 
the stellar emission (where the visibility V* = 1 as the predicted 
angular sizes of the stars are 0.63mas for HD69830 and 0.73mas 
for j] Corvi based on typical radii for their spectral types and 
their distances determined by parallax) and the component from 
the disc (where the Vdi sc will depend on the resolution of the disc 
flux). The total visibility of a star + disc source will then be given 
by 

Vtot = |^V* + §^Vd lsc , (3) 
r tot r to t 

(where ,F+ is the flux of the star, etc.). As only the disc com- 
ponents of our targets have the potential to be resolved on the 
baselines used in this study the equation above predicts a visibil- 
ity function comprising the sum of two terms. The first of these 



3.1. HD69830 

HD69830 was observed on all four baseline combinations used 
for this study. A sketch of the geometry of the baselines and 
relative lengths is shown in Figure|4](top left). 

The calibrated correlated flux for all four baseline observa- 
tions of this source are shown in Figure [4] (top right). The er- 
ror bar indicates an average error measured across all four base- 
lines. These errors are determined from two sources: splitting the 
fringe tracking files into 5 equal subsets and performing the re- 
duction on each subset in turn to determine the variation in fringe 
signal across the time of the observation (statistical error); and 
the variation between correlated flux measurements as calibrated 
by different standard star observations where available. The two 
sources were added in quadrature. The error is dominated by cal- 
ibration error which is typically ~10% of the correlated flux (see 
also 2.4). 

The correlated fluxes show variation across the different 
baselines, as can be seen in the levels of flux measured at 8 pm. 
When correlated fluxes are different on different baselines there 
is the possibility that these differences represent real differences 
in the resolution of the source. However, at 8 pm the disc emis- 
sion is expected to be very low (4 mJy at 8pm), and so ~ 100% of 
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Fig. 4. The results of MIDI observations of HD69830. Top left: The baselines observed in this program. Top right: The calibrated 
correlated fluxes found in the MIDI observations of HD69830. The thick line represents the total IRS spectra, which represents the 
maximum flux we should see. Note that although the absolute values are varying by ~ 10%, all runs show a fairly flat spectrum, with 
no evidence of the silicate features seen in the IRS spectra. Bottom left: The visibility of HD69830. The measured visibilities on all 
four baselines are calculated using the IRS spectra and the measured correlated fluxes. The prediction for what the visibility curve 
should look like if the disc flux is totally resolved (V res , equation|4| is also shown by a thick solid line. Bottom right: The visibility 
of HD69830 scaled to an average of 1 between 8-9 yt/m. Notice the shape across all baselines is similar to the V res prediction. 



the emission at 8/mi should come from the stellar photosphere. 
As the star is expected to be point-like and completely unre- 
solved, the correlated flux should equal the total flux (marked 
by a thick solid line) at this wavelength. The differences in these 
absolute values are at the expected level of variation in absolute 
correlated flux from changing conditions between science and 
standard star observations (see section 2.4). 

What is clear from Figure [4] (top right) is that for all base- 
lines the correlated flux is flatter than the total flux measured in 
the IRS spectrum; there is no evidence of the silicate feature seen 
between 10-1 1.5 /urn in FigureQ] As shown by equation|4j if the 
disc flux is totally resolved (V res ) the correlated flux should con- 
sist solely of the stellar photospheric emission - with no evidence 
of the silicate features seen in the excess emission spectrum. The 
indication that the disc flux has been at least partially resolved 
on all baselines is made clearer by examination of the visibility 
function shown in Figure [4] (bottom left). 

The absolute values of visibility across all baselines are an 
average of 1.04±0.16. That this variation is a reflection of abso- 
lute uncertainty and not evidence of different levels of resolution 
is confirmed by the visibility at 8 yum, which as discussed above 
should be 1 as the flux here is completely dominated by the stel- 
lar photosphere. Uncertainty at this level is to be expected, given 
typical uncertainties in calibrated visibilities of 0.10 (see Figure 
[3]) and the relative faintness of the target (total flux of ~1 Jy 



across the MIDI spectral range, see Figure [TJ. Beyond this un- 
certainty in absolute values, however, the visibilities measured 
on all four baselines can be seen to dip exactly in the regions we 
would expect a dip if the disc flux was resolved. This dip is made 
clearer if we scale the visibility function to 1 averaged over the 
short wavelength region (8-9/mi), as would be expected even in 
the resolved disc case (V KS ) because of the lack of excess in this 
region. 

The scaled visibility functions measured on the four base- 
lines are shown in Figure [4] bottom right. The dip in visibility 
between 10-1 1.5 fim is now very obvious. To confirm the signif- 
icance of this dip, the unsealed average visibility within the 14 
spectral channels between 8-9 yum was compared to that in the 
29 spectral channels between 10-11.5/vm (i.e., around the low- 
est predicted visibility when Vdi sc = 0). Table [3] presents these 
weighted means (the weights were determined from the statis- 
tical variation within the 5 sub-integrations available for each 
spectral channel) together with their standard errors. Calibration 
uncertainty is not included in these values. The visibility ratios 
presented in the table are the ratio of these weighted means, 
where we have included a contribution from calibration errors in 
quadrature with the statistical errors in the 2nd and 3rd columns 
when evaluating the errors in the final column. These calibra- 
tion errors arise from two sources. First, as we have used the 
IRS spectrum to determine the visibilities, we included an error 
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Baseline 




Visibility 




n£imc 


8-9/mi 


10-11 Sum 


Ratio 


A 


0.944 ± 0.003 


0.879 ± 0.005 


0.931 ± 0.037 


B 


1.003 ±0.003 


0.928 ± 0.004 


0.926 ± 0.027 


C 


1.093 ±0.003 


1.018 ±0.004 


0.931 ±0.028 


D 


0.839 ± 0.006 


0.803 ± 0.005 


0.958 ± 0.039 


v res 


0.990 


0.860 


0.869 



Table 3. Measuring the significance of the dip in visibility of 
HD69830 at 10-11.5 /mi. The values of visibility in each wave- 
length range are the weighted means taken over all spectral chan- 
nels in that range, and uncertainties quoted are the standard er- 
rors on the weighted mean from statistical uncertainty only (un- 
certainty in the absolute visibility is higher due to calibration un- 
certainties). Final error on ratio includes error from the calibrator 
ratios discussed in section 2.4 and a 1.3% error on the spectral 
slope of IRS observations (Beich man et al.l l2006). as well as the 
statistical errors from each wavelength range. The V res predicted 
visibility arises from assuming that Vdi sc = (see equation|4]i. 



of 1.3% typical of errors on the IRS spectral slo pe as measured 
for so urces with no observed excess emission feeichma n et alj 
2006). Second, we included the error on the ratio from calibra- 
tion as measured in the standard-standard calibrated visibilities 
discussed in section 2.4. We used the standard deviation of the 
visibility ratio (3.7%, section 2.4), but for runs B and C where 
2 standard star observations were available for calibration, this 
was reduced by a factor of V2. 

For all four baselines used to observe HD69830 the visibility 
ratio suggests a partial resolution of the disc - the ratio is higher 
than the 0.869 expected for V res but definitely lower than unity. 
On baselines B and C this difference from unity is significant at 
a level of ~ 3<x. When averaged over all 4 baselines the weighted 
mean of the visibility ratio is 0.934 ±0.0150 Taken as a whole, 
these data suggest that the observed visibility ratio is different 
from 1 at a level of 4<x. We note in passing that the actual scat- 
ter between the four measurements of the visibility ratio is very 
small, suggesting that we have been conservative in associating 
an error of 0.015 to our final "best-estimate". 

We do note, however, that the lower right panel of Figure 
[4] also shows the visibility function of HD69830 rising beyond 
~ 12//m to above 1 in the scaled case, meaning the ratios between 
the weighted mean visibilities between 1 1.5-1 3//m and between 
8-9^m are on average greater than 1 (baseline A 1.045, B 0.991, 
C 1.043 and D 1.098). For baseline D this ratio is larger than 
would be expected from any calibration error, as the error on 
the ratios over the same wavelength ranges from the calibrator- 
calibrator visibilities (section 2.4) is a maximum of 0.053. The 
V res prediction shows that the visibility function in the case that 
the disc is completely resolved would be expected to rise, but 
visibilities should not ever be greater than 1 (see equation 0). 
This is unlikely to be due to the water ice discussed in sec- 
tion 2.3, as if some of the excess emission measured in the IRS 
spectrum did fall outside of the MIDI beam we would be over- 
estimating 7 tot in our calibration, and this would lead to a falsely 
low visibility. We consider that the lower signal-to-noise in the 
MIDI fringe measurements on this target at longer wavelengths 
(signal-to-noise is an average of ~70 between 8-9yum, 37 be- 



1 Note that these calculations assume that the visibility on all base- 
lines is the same. Although this is not generally expected to be the case 
it is true if the disc is completely resolved or unresolved and so these 
calculations are used to determine the significance with which we can 
say that the disc is completely resolved or unresolved. 



tween 10-11.5//m and 15 beyond 12yum) is the cause of these 
high visibilities which causes a residual slope. Whether this af- 
fects the 10-1 1 .5yum region is not clear but for this reason we do 
not consider that the results provide strong evidence for a par- 
tially resolved rather than completely resolved disc. The limits 
we can place on the morphology of the HD69830 debris disc 
with the results presented in this section are discussed in section 
5.1. 

3.2. r] Corvi 

r\ Corvi was observed on baselines A and B and on baseline D 
twice (see Figure [5] top left panel). The second observation on 
baseline D was carried out with an increased chop frequency and 
integration time which led to greatly improved intensity calibra- 
tion data. 

The calibrated correlated fluxes for all baselines are shown 
in Figure|5] top right, where we have included calibration errors 
and the variation across sub-exposures added in quadrature. The 
error bar shown is indicative of the average error observed over 
all baselines. Overpl otted on these figu res is the total IRS pho- 
tometry presented in lChen et al. (2006), as also shown in Figure 
ffl 

The correlated fluxes are very similar in shape across all 
four observations. The absolute values vary by the typical ~10% 
we expect for absolute variation in correlated flux (see section 
2.4), but all show a spectrum consistent with photospheric emis- 
sion in the Rayleigh-Jeans regime. In common with the result of 
HD69830, there is no evidence of the spectral feature seen in the 
IRS total spectrum in any of the correlated flux measurements. 
This indicates that for this source the disc emission also appears 
to be partially resolved. Again, the result is made clearer by ex- 
amination of the visibility function. 

The calibrated visibilities for all observations of r\ Corvi are 
shown in Figure [5] (bottom left), along with the visibility pre- 
dicted for a completely resolved disc component marked in a 
thick solid line (V res , see equation 0). The dip in visibility at 
~1 1.5 yum is seen very clearly on all observed baselines. Again, 
as expected, the variations between different baselines in terms 
of absolute values are an average of ~0.1, with a maximum dif- 
ference of ~0.2. The visibilities scaled to the V res prediction be- 
tween 8-9//m are shown in Figure|5](bottom right). Although we 
do not know what the visibility should be at 8-9//m for r\ Corvi 
(as we do for HD69830) as there is significant disc flux in this 
range, this scaling allows the shape of the visibility function to 
be seen more clearly and allows easy comparison with the V res 
model. 

In Table 2] the weighted mean visibilities in the short wave- 
length region (8-9 //m) and the region in which the silicate ex- 
cess feature has the strongest contribution to the total flux (10- 
11.5 fim) are listed. Standard errors and errors on the ratio in- 
cluding calibration uncertainty are calculated as described for 
HD69830 in section 3.1. As is clear from the table the differ- 
ence seen in the measured visibilities at shorter wavelengths and 
midway through the MIDI range (10-1 1.5yum) is significant. Not 
only this, but the size of the change is similar in all cases to 
that predicted for the case of completely resolved excess emis- 
sion, as can be seen in the scaled visibility function plots (Figure 
[5] bottom right). Results on individual baselines are incompati- 
ble with an unresolved disc (ratio = 1) at a level of at least 4 <x 
(see Table|4|. Taking all the results together (see footnote[Tji the 
weighted mean dip in visibility is 0.880 ± 0.013, incompatible 
with an unresolved disc at the 9<x level. The individual measured 
ratios and average result are close to the prediction for a com- 
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Fig. 5. The results of MIDI observations of 77 Corvi. Top left: The baselines observed in this program. Top right: The calibrated 
correlated fluxes found in the MIDI observations of 77 Corvi. The thick line represents the total IRS spectra, which represents the 
maximum flux we should see. Note that although the absolute values are varying by ~ 10%, all runs show a fairly flat spectrum, with 
no evidence of the silicate features seen in the IRS spectra. Bottom left: The visibility of 77 Corvi. The measured visibilities on all 
four baselines are calculated using the IRS spectra and the measured correlated fluxes. The prediction for what the visibility curve 
should look like if the disc flux is totally resolved (V res , equation|4| is also shown by a thick solid line. Bottom right: The visibility 
of 77 Corvi scaled to the average visibility from the V re .s prediction between 8-9 fim. This scaling shows the four observations all 
show strong evidence of being similar to the V les prediction. 



pletely resolved disc (V les ). Note that we do not see the same ev- 
idence for a significant residual slope in the visibility function as 
was the case for HD69830. The observations of 77 Corvi benefit 
from a higher signal-to-noise even at longer wavelengths (aver- 
age signal-to-noise in the correlated flux measurements is 128 at 
8-9/um, 75 at 10-1 1.5/jm and 31 at wavelengths beyond 12/jm). 
Though the levels of absolute calibrated flux and thus visibil- 
ity are not constrained to better than 10%, (see section 2.4) the 
dip and its appearance at a significant level across all baselines 
is consistent with the disc emission being completely resolved 
on all observed baselines. The constraints we can place on the 
morphology of this source based on these results are discussed 
in section 5.2. 



Baseline 




Visibility 




name 


8-9/im 


10-11.5/jm 


Ratio 


A 


1.016 ± 0.005 


0.892 ± 0.004 


0.877 ± 0.026 


B 


0.910 ±0.003 


0.818 ± 0.003 


0.899 ± 0.027 


Dl 


0.841 ± 0.003 


0.739 ± 0.003 


0.879 ± 0.026 


D2 


0.922 ± 0.003 


0.801 ± 0.004 


0.868 ± 0.026 


v res 


0.915 


0.784 


0.857 



Table 4. Measuring the significance of the dip in visibility of 
77 Corvi at 10-11.5 //m. A description of the calculation of the 
visiblities and errors is given in the caption to Table [3] 



4. VISIR imaging of HD69830 

Following t h e 8m limits on the 77 Corvi emission presented in 
ISmith et all d2008l) . we have new data obtained on VISIR at 
the VLT under proposal 079.C-0259 on 6 th April 2007 on the 
emission around HD69830. The data were all taken in filter 
Q2 {A c = 18.72/um, AA = 0.88/jm) with pixel scale 07075. 
The observations were performed using a chop and nod throw 
of 8". Chopping was performed in a North-South direction 
and nodding was performed in the perpendicular direction. The 



observations consisted of two integrations of 2100s each on 
HD69830, with standard star observations before and after each 
science target integration (125s per integration). Frequent stan- 
dard star observation was used to determine how the PSF was 
varying through the observations. The standard star observa- 
tions were also used for photometric calibration. The standard 
star, HD61935, was cho sen from a list of m id-infrared spectro- 
photometric standards dCohen et alj 1 19991) and was also used 
in the MIDI observations (Table [T]). The observations are sum- 
marised in Table [5] 
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Fig. 6. The VISIR images of HD69830 and standard star observed in filter Q2. Left: The final image of standard star (HD61935) 
after co-addition of 4 sub-integrations. This image acts as a PSF reference. Middle: The final image of HD69830. Contours on this 
and the standard star are at levels of 25, 50 and 75% of the peaks. No evidence of emission extended beyond the PSF reference 
is seen in the image of HD69830. The colour scale on these images goes from to the peak. Right: The image of HD69830 after 
subtraction of the PSF reference (standard image) scaled to the peak. The colour scale has a minimum of and maximum of 3<x 
per pixel. Contours mark any region of >2sigma significance. The stripe to the South East of the image is a detector artifact. No 
significant emission is seen after subtraction of the PSF which would indicate resolved emission (testing over a range of regions of 
different shapes which could detect a range of different disc geometries gives a maximum detection of <2<x in this image). 



Target 


Target type 


Integration time (s) 


HD61935 


Cal 


125 


HD69830 


Sci 


2100 


HD61935 


Cal 


125 


HD61935 


Cal 


125 


HD69830 


Sci 


2100 


HD61935 


Cal 


125 



Table 5. The observations of HD69830 and standard star with 
VISIR on the VLT. Observations were performed under proposal 
079.C-0259. All observations were performed using filter Q2 
and the small field detector (pixel scale (X'075). 



The d ata were reduced u sing custom routines described in 
detail in ISmith et aU d2008l) . In summary, data reduction in- 
volved determination of a gain map using the mean values of 
each frame to determine pixel responsivity (masking off pixels 
on which emission from the source could fall, equivalent to a 
sky flat). In addition a dc-offset was determined by calculating 
the mean pixel values in columns and rows (excluding pixels 
on which source emission was detected) and this was subtracted 
from the final image to ensure a flat background. Pixels showing 
high or low gain, or those which showed great variation through- 
out the observation, were masked off. The chop-and-nod pattern 
adopted results in four images of the star falling on the detec- 
tor. These images were co-added after determining the centre 
of each image by 2-dimensional Gaussian fitting and then the 
sub-integrations were added together to give a final image of 
HD69830 and the standard star which would act as a PSF refer- 
ence. 

The photometry was performed using circular apertures of 
1 "radius centred on the peak of the emission. Statistical back- 
ground noise was calculated in annuli centred on the peak of 
emission with inner radius 2", outer radius 4". Calibration un- 
certainty was estimated from the variation in calibration levels 
from the four sub-integrations on the standard star, and was 
found to be 9%. HD69830 was detected in the final co-added 
image with a signal-to-noise of 32, with the total photometry in- 
cluding calibration uncertainty of 377 + 46 mJy. The IRS spec- 
trum of this source has a flux of 365 mJy at 18.72/vm (253 mJy 



expected from the photosphere). No colour-correction was ap- 
plied for this narrow-band filter. In addition, no correction for 
airmass was employed, as the standard star was observed at a 
similar airmass to HD69830 and no evidence was found that the 
calibration factors from the sub-integrations on the standard star 
were correlated with airmass. 

The final images of the standard star HD61935, used as 
a PSF reference, and of HD69830 are shown in Figure [6] 
Overplotted on the images are contours at 25, 50 and 75% of 
the image peak (brightness peaks were 7980 mJy/arcsec 2 for 
the standard star observation, 1033 mJy/arcsec 2 for HD69830). 
Comparison of the contour plots indicates that there is no sign of 
extended emission around HD69830 that could be indicative of a 
resolved disc. This is confirmed by examination of the residuals 
image (Figure|6]right). This is the HD69830 image after subtrac- 
tion of the PSF reference image (HD61935, Figure|6]left) scaled 
to the peak of emission. The residual image shows no evidence 
for emission beyond the PSF. Further tests comparing the pro- 
files of the observations, Gaussian and Moffat profile fits to the 
images and comparisons of FWHM measurements all confirmed 
this result. Different test regions (annuli of different dimensions 
for face-on disc geometries, rectangular boxes for edge-on discs, 
and intermediate shapes for discs with intermediate inclinations 
to the line of sight) were also examined for any evidence of sig- 
nificant emission in the residuals image. No such emission was 
found (most significant positive flux of just less than 2cr signif- 
icance was found in a region including a detector artifact - see 
Figure |6). 

We followed the procedure of ISmith et al.l d2008l) to place 
limits on the size of the disc after no extension was dete cted. 
This procedure is described in detail in lSmith et al.l ((2008), and 
involves convolving the PSF model with different disc + star 
models to determine what combinations of disc parameters (ge- 
ometry and emission level) would have resulted in a detection 
of extended emission. The limits depend strongly on the sen- 
sitivity of the observation and the variation of the PSF, which 
was determined in this case from examination of the four sub- 
integrations on the standard star. The PSF was found to have a 
FWHM of 0752±0.04, a variation of 7%. The individual inte- 
grations gave four separate PSF models which were used in the 
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Fig. 7. The limits we can place on the geometry of the HD69830 
following no detection of extended emission in the VISIR Q 
band imaging. Discs in the space above the lines should have 
been detected in our imaging, and those lying in the shaded re- 
gion would have remained undetected regardless of geometry. 
Different lines represent different disc geometries as given in the 
legend. Narrow discs have a disc width of 0.2r, wide discs a 
width of 2.0r, where r is the radius of the mid-point of the disc. 
The expected disc flux and location from the SED is marked by 
an asterisk, with error bars representing the uncertainty in R \ 
from errors in the IRS spectrum. 



limits testing, with variations between the results being the limit- 
ing factor on disc detection at small disc radii. The results of this 
testing are shown in Figure[7] For a face-on disc assuming a nar- 
row ring-like structure (width Q.2r where r is the disc radius) the 
extension testing places a limit of <0719 ± Of! 01 on the disc ra- 
dius (<2.4 AU at a distance of 12.6pc). This limit is determined 
by the size and variation of the PSF (limiting line at small disc 
radii is near vertical as it is not strongly dependent on the disc 
flux). A PSF variation of 7% is fairly typical for Q band observa- 
tions taken under g ood conditions based on previous observing 
experience (see e.g. lSmith et a l. 2008, 2009). Only observations 
taken under optimal (very stable seeing) conditions would allow 
us to probe closer to the star, and then only by perhaps one or 
at best two tenths of an arcsecond (~0. 1-0.25AU). This limit as- 
sumes disc emission is 112±15mJy (3<x limit using errors from 
IRS photometry), as measured in the IRS spectrum assuming a 
Kurucz model profile for the photosphere and in agreement with 
our photometry. Limits on broader discs arise from the sensitiv- 
ity of the observations, and could be improved with increased 
observation time. With the current limits, we can rule out any 
disc brighter than 33mJy at a radius of Of! 5 (6.3AU) with some 
dependence on the assumed disc geometry (Figure|7|. 

5. Discussion 

We now consider the implications for the geometries of the 
HD69830 and rj Corvi discs in the mid-infrared following the 
MIDI results and the constraints from 8m imaging. 

5.1. HD69830 

The MIDI correlated fluxes measured for HD69830 and the vis- 
ibility functions determined from these are consistent with the 
detection of photospheric emission only. The slopes are consis- 
tent with a Rayleigh-Jeans slope, and there is no evidence of the 



silicate feature observed in the IRS spectrum. We thus conclude 
that the excess emission does not strongly contribute to the cor- 
related flux and therefore this has been at least partially resolved 
on the observed baselines. 

Evidence for resolved emission around £ Aql in the near- 
infrare d was interpreted as evidence for a low-mass companion 
object dAbsil et alj|20 08). We do not believe this is the explana- 
tion for the resolved emission around HD69830, since the emis- 
sion spectrum is consis tent with dust emission (Beichma n et alj 
l2005tlLisse et a l. 2007), and the observed visibility ratio is the 
same (within the errors) on all observed baselines. The same 
resolution on near perpendicular baselines such as A and B/D 
(see Table [2] and Figure [4] top left) would not be expected if the 
emission arose from a point-like source. Furthermore extensive 
radial ve locity observatio ns have revealed only low-mass com- 
panions dLovis et alJl2006l) orbiting within 1AU, believed to be 
Neptune-mass (although the masses could be larger if the sys- 
tem is observed face-on). Such planets cannot be the source of 
the observed change in visibility with wavelength due to the low 
expected stanplanet mid-infrared fl ux contrast of < 10~ 3 (see e.g 
iBurro ws 20051 iBurrows et alj|2006l) . 

The observed visibilities are consistent with a partial resolu- 
tion of dust emission around HD69830. Assuming that the stel- 
lar component has a visibility of 1, and that the values of and 
Fdi S c are accur ately given by a Kuru cz model profile and the IRS 
spectra from Beichman et alj d2005l) . the visibility ratios listed 
in Table [3] can be used to estimate the disc location given certain 
assumptions about the morphology of the emission^ For exam- 
ple, assuming that the disc emission has a Gaussian distribution 
gives a best fit FWHM of lOmas. As many debris discs resolved 
to date have been found to have ring-like structure (see, e.g. 
iGreaves et al]2005tlKalas et alj2005tls"chneider et al.l2009h . we 
also modelled the observed visibilities assuming the disc emis- 
sion to originate in symmetric, face-on rings with varying radii 
and widths of dr/r = 0.2, 1 or 2 where r is the radius of the 
middle of the ring (see appendix A). Figure [8] shows the results 
of a ^f 2 -test of the goodness-of-fit of each model suggesting a 
best-fit disc radius of 10-25mas, depending on the width of the 
ring. Both tested morphologies predict a best-fit radius of 0.1- 
0.3AU at 12.6pc which would place the dust in the middle of the 
planetary system. A stability analysis bv lLovis et alj {2006) used 
Monte-Carlo modelling to determine locations where massless 
particles could have long-term stable orbits in the presence of 
the 3 Neptune-mass planets, finding 2 long-term stable zones at 
0.3-0.5 AU and >0.8AU. Thus the visibilities suggest a best-fit 
location for the dust compatible with the inner stable zone. 

However, uncertainties in the visibility ratios mean that the 
observations are also consistent within the uncertainties with a 
bigger radial location (see Figure [8]), particularly if the possible 
slope on the visibilities is taken into account (see section 3.1). 
Should the high levels of visibility seen at long wavelengths be 
evidence of a residual slope, then the disc emission may be com- 
pletely resolved on these baselines (in which case the visibil- 
ity would match the V res prediction). For this to be the case we 
would require a Gaussian of FWHM > 22mas (0.28AU) on an 
average baseline of 83m, or of at least 40mas (0.5AU) to be com- 
pletely resolved on baseline B ( 44m). These sizes a re consistent 
both with the inner stable zone o flLovis et alj d2006l) . and a larger 
radial location outside the orbits of the planets. A larger radial 
location is favoured from analysis of the IRS spectrum which 



2 A detailed description of how the models were compared to the 
observed visibility ratios and how errors were propagated through the 
analysis is given in Appendix A. 
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Fig. 8. The limits on the HD69830 disc from the MIDI observations. Left: The confidence limits found comparing the observed 
visibility ratios on all 4 baselines with those predicted by face-on ring-like models for the disc morphology (of radius r and width 
dr). Lower values of x 2 represent a better fit to the data, and the horizontal dotted lines denote the level at which the models are 
consistent with the observations at the 1, 2 and 3<x levels as determined by the percentage points of the^ 2 distribution function with 
4 degrees of freedom. The limits depend on the assumed width of the ring as indicated in the legend. Right: Predicted visibility 
ratios for narrow (dr/r = 0.2) ring-like models of radii 5, 10, 25 and 80mas (solid lines). The ratio expected if the disc emission is 
completely resolved (i.e., the V res prediction) is shown by a dashed line. Observed visibility ratios are shown with asterisks and lcr 
error bars (TableO. The modest level of significance with which different models fit the observations (seen in the left figure) is due 
to the complex visibility function for ring models, combined with sparse baseline coverage and large uncertainties. See Appendix 
A for details of modelling and error calculations. 



suggests the emission originate s in disrupted P o r D-type aster- 
oids at a distance of 0.9-1 . 1 AU dLisse et al.ll2007l) . If there is not 
a residual slope on the data, and the results are truly indicative 
of a partial disc resolution, the dust would then be constrained 
to lie within the orbits of the planets challenging current SED 
models. 

Regardless of the residual slope, the MIDI observations are 
inconsistent with unresolved emission, and so limits can be 
placed on the minimum radius of the dust emission. Using ring- 
like models we find a 3<x limit on the radius of >4mas (0.05AU), 
or from Gaussian models we find a 3cr limit on the FWHM of 
>4mas. The VISIR data (section 4) provides a complimentary 
upper limit on the disc radius of <2.4AU. 

Adopting the radial size of 1 AU from the SED modelling, we 
compared the observations to model predictions for 1AU rings 
inclined both edge-on, face-on and 45° to the line-of-sight and at 
different position angles to assess whether the observations show 
evidence for an inclined disc structure; measurement of the disc 
inclination would remove ambiguity in the planet masses assum- 
ing these were orbiting in the same plane. A narrow edge-on disc 
lying at a position angle of ~165° provides the best fit to the four 
observed visibility ratios (y 2 minimum of 1.3<x, Figure|9j other 
disc widths provide worse fits to the data, albeit still within ob- 
servational uncertainties). Inclined discs at other position angles 
away from ~60° also provide better fits to the data than a face- 
on model, although a narrow face-on disc at 80mas is consistent 
with the observed visibilities at a level of 3<x. Thus with the cur- 
rent data we cannot place significant limits on the disc inclina- 
tion. The lack of significant differences in the visibility ratios on 
different baselines also means there is no evidence for clumpy 
and asymmetric structure, although the uncertainties in the data 
mean that such structure cannot be ruled out. 

Taking the stab ility analysis of lLovis e t al. (2006) and the 
SED modelling o flLisse et alid2007l) together with the MIDI and 
VISIR results presented here, we consider the emission from this 
system to be most likely smooth and symmetric at ~1AU radius. 
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Fig. 9. Comparison of inclined disc models to the observed 
visibility ratios of HD69830 at a fixed disc radius of 80mas. 
Inclinations are measured from face-on = 0°. 



5.2. rj Corvi 

The MIDI correlated fluxes are consistent with Rayleigh- Jeans 
emission across the MIDI range. We see no evidence of the mid- 
infrared excess emission dominated by silicate emission in the 
correlated fluxes and conclude that the data are consistent with 
all of the excess being resolved on the observed baselines. As 
for HD69830, we rule out a companion object as the source of 
the resolved emissio n due to the silica te emission feature seen 
in the IRS spectrum dChen et al.| [2006) and the lack of a known 
binary companion. Limits from radial velocity searches rule out 
a planet of mass >0.4Mj up with a period of 3 days or less, and a 
planet of mass >2.1Mj llp with a period of 100 days or less (cor- 
respo nding to <0.46AU for M+ = 1.5M ), at the 3<x confidence 
level dLa grange et al .1120091) . 

The observed visibilities are consistent with at least par- 
tial, perhaps complete, resolution of the dust emission. The ob- 
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Fig. 10. The limits on the rj Corvi warm dust from the MIDI observations. Left: The confidence limits found comparing the observed 
visibility ratios on all 4 baselines with those predicted by face-on ring-like models for the disc morphology (of radius r and width 
dr). Lower values of^- 2 represent a better fit to the data, and the horizontal dotted lines denote the level at which the models are 
consistent with the observations at the 1, 2 and 3<x levels as determined by the percentage points of the^ 2 distribution function with 
4 degrees of freedom. The limits depend on the assumed width of the ring as indicated in the legend. Right: Predicted visibility 
ratios for narrow (dr/r = 0.2) ring-like models of radii 5, 10, 25 and 90mas (solid lines). The ratio expected if the disc emission is 
completely resolved (i.e., the V res prediction) is shown by a dashed line. Observed visibility ratios are shown with asterisks and lcr 
error bars (Table|4]i. The modest level of significance with which different models fit the observations (seen in the left figure) is due 
to the complex visibility function for ring models, combined with sparse baseline coverage and large uncertainties. See Appendix 
A for details of modelling and error calculations. 



served visibility ratios were modelled in a similar manner to 
HD69830. Assuming that the disc emission has a Gaussian dis- 
tribution results in a best fit FWHM of 27mas. The compari- 
son with ring-like models is shown in Figure \TU[ from which 
best-fit disc radii are 20-25mas for wider rings, or >50mas for 
narrow face-on rings (>0.91AU). A size of 20-50mas translates 
to a radius of 0.4-0.9AU at 18.2pc. Since the MIDI observa- 
tions are inconsistent with unresolved emission, they can be used 
to derive a minimum size for the emitting region. Assuming a 
Gaussian distribution for the disc emission limits the FWHM to 
>14mas (>0.25AU; 3cr limit from x 2 goodness-of-fit over all 4 
baselines), whereas the ring-like models constrains the radius to 
>9mas (0.16AU, 3cr limit). 

There remains some uncertainty from SED modelling on the 
location of the mid-infrared excess, with two models proposed 
that suggest that either the hot dust lies in a si ngle ring at 1 ,7AU 
or in two radial lo cations at 1.3 and 12AU dChen et alj |2006). 
ISmith et all d2008l) ruled out the 12AU component at the 2.6<x 
level, since this was not resolved in VISIR 18.7/im imaging, so 
that a single ring at 1 .7AU remains the best estimate of the radius 
of the hot dust from the SED, with unresolved 8m imaging con- 
straining the location of a single radius component to < 3AU. 
Thus, like HD69830, the best-fit radius is smaller than, but not 
inconsi stent with that pr edicted from SED modelling (90mas, 
1.7AU,lSmidieLalj2008); the narrow ring models are most con- 
sistent with the SED models, however the errors on the data and 
the complex visibility functions expected for ring-like structure 
(FigureQI)] right) mean that a large range of ring models provide 
a good fit to the data within the errors. In contrast to HD69830 
no planets are known to be orbiting interior to th e dust, although 
limits are not so stringent dLagrange et alj|2009l) . 

Adopting a radius of 1.7AU for the dust from SED mod- 
elling, we compared the observations to model predictions for 
1.7AU rings inclined both edge-on, face-on and at 45° and 70° 
to the line-of-sight and at different position angles to attempt 
to constrain the inclination and position angle of such a disc. 
The best fitting disc geometry was found to be a broad edge- 



on disc at a position angle of ~120° (Figure ITTb. For narrow 
discs all geometries are within l<x of the best fitting model, and 
for broad discs only edge-on discs with a position angle ~15° 
are outwith 1 cr of the best fit. Thus other disc geometries can- 
not be ruled out with the current data. However, it is notable 
that this best-fit position angle coincides with that of a cold disc 
component (at 130±10°) imaged in the sub-mm and inferred 
to be in a 150+20AU r ing inclined by 45 ±25° to our line-of- 
sight dWvatt et alj |2005). A common position angle would sug- 
gest that both dust components, over two orders of magnitude in 
radius, originate from planetesimals orbiting in the same plane. 

There is no evidence from this data for asymmetric of 
clumpy structure, since the visibility functions are consistent 
(within the errors) with the same level of resolution on all base- 
lines (Figure |5). This level is consistent with the excess being 
completely resolved on all baselines (i.e., it coincides with the 
Vres prediction). We note that this does not mean that the emis- 
sion comes from a region which appears the same size on all 
observed baselines, but that on each baseline the emission cov- 
ers a large enough region to have been completely resolved. 

Taking the MIDI results and c ombining them with the results 
from 8m imag ing in ISmith et al.l d2008l) and the IRS spectrum 
presented in IChen et aTl d2006l) . we consider the emission from 
this system to be most likely smooth and symmetric at 1 .7AU ra- 
dius, with tentative evidence for extension along a position angle 
aligned with the cool outer disc. 

5.3. Transient sources of emission 

For both sources we infer the dust to lie at 1-1 .7AU and so con- 
sider that it must have a transient origin dWvatt et al.ll2007l) . One 
possible source for the transient emission is a recent massive 
collision within a coincid ent planetesimal belt (a possibility pro- 
posed for BD+20307, see lZuckerman et al.l2 008, and references 
therein). The emission from such an event might be expected to 
be clumpy or asymmetric as the dust would start off concentrated 
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Fig. 11. Comparison of inclined disc models to the observed 
visibility ratios of r\ Corvi at a fixed disc radius of 90mas. 
Inclinations are measured from face-on = 0°. 



in a small spatial region (rather than dispersed in a ring like an 
asteroid belt for example). There is no evidence for such a struc- 
ture from the MIDI observations. Furthermore this scenario was 
considered as possible but unlikely due to the expected low fre- 
quency of sufficiently massive collisio ns and the lifetime of re- 
sulting dust particles from such events ( Wyatt et al. 2007). Thus 
it seems more likely that the parent bodies of the dust seen at 
~1AU originated in a planetesimal belt at larger distances where 
collisional timescales are much longer. The mechanism trans- 
porting the dust into 1AU remains unknown, though one possi- 
bility is a dynamical instability in a period analogo us to the Late 
Heavy Bombardment (L HB) in the Solar System (Gome set alj 
120051: iBooth et"aT]|2009t) . the cause of which could be a recent 
stellar flyby or the migration of giant planets in the sy stem. 

The known cool dust population toward r\ Corvi dWvatt et alj 

120051) could represent the outer parent planetesimal belt sus- 
pected to be feeding the hot dust, a suggestion supported by the 
tentative evidence from the MIDI observations for a position an- 
gle for the hot component (120°) that is similar to that of the 
cool dust (130±10°). However, the errors and sparse baseline 
coverage of the current data do not allow significant limits to 
be placed at this stage on the hot disc's inclination and position 
angle. In contrast no cold dust population has yet been detected 
around HD69830. Indeed the MIP S photometry limits th e ex- 
cess emission at 70//m to 1+3 mJy (Beichman et al. 2005), with 
SCUBA photomet r y limiting excess flux at 850/vm to <7mJy 
dSheret et al.ll2004t [Matthews et"ai]|2007l) . indicating that there 
is little if any cold dust in the system. One possible way around 
th e conclusion that th e dust has to be transient was suggested 
bv lPavne et alj d2009h . who considered the dynamical evolution 
of the HD69830 pla netesimal disc du ring the growth and migra- 
tion its 3 planets (Alibert et al. 2006). In their model significant 
planetesimal mass remained outside 1AU on high eccentricity 
orbits, and they suggested that the increased collisional lifetime 
of highly eccentric orbits could allow significant mass to persist 
even to 2Gyr, a possibility which is currently being investigated. 

6. Conclusions 

In this paper we present MIDI observations of HD69830 and r\ 
Corvi, which for the first time resolve the hot dust emission in 
these systems, and setting lower limits on the radial location of 
the emission of >0.05AU and >0.16AU, respectively. We also 
present new VISIR observations of HD69830 which show no 



evidence for resolved emission, thus constraining the radial lo- 
cation of the emission to <2.4AU. Similar 8m observ ations of 
r\ Corvi were already presented in ISmith et alj (2008) showing 
this emission originates <3AU. These mid-infrared observations 
thus place the dust squarely at a location consistent with that 
predicted from SED modelling at 1.0 and 1.7AU for HD69830 
and r\ Corvi respectively. At such a location the emission to both 
sources is expected to be transient. 

There is no evidence from the MIDI data that the emis- 
sion is anything but smooth and symmetric. However, assum- 
ing a 1.7AU radius for the emission toward r\ Corvi the data 
marginally favour a more edge-on disc lying at a position angle 
consistent with that of the known 150AU cool disc at 130±10°. 
Such a correlation would further support the suggestion that the 
cool outer disc is feeding the inner disc. The MIDI data also 
suggest that the HD69830 emission is partially, rather than com- 
pletely, resolved, placing it at 0.1-0. 3 AU. Observational uncer- 
tainties, including a possible residual slope, mean that this inter- 
pretation is not favoured at present given the SED constraints, 
but if confirmed would present a challenge to the current mod- 
elling of the IRS spectra of this source. 

These observations highlight the viability of exploration of 
dust in terrestrial planet regions using mid-infrared interferom- 
etry. Future MIDI observations of HD69830 and r\ Corvi that 
extend the baseline coverage would allow more stringent con- 
straints to be placed on the geometry of these debris discs, in- 
cluding measuring the disc inclination for HD69830 (thus con- 
straining the masses of planet found by radial velocity studies), 
the disc position angle for r\ Corvi (for comparison with that of 
the outer cool disc), and searching for the clumpy or asymmetric 
structure that might be expected from a recent massive collision. 
Such observations thus have significant potential to constrain the 
different theories for the origin of these unusual hot dust popu- 
lations. 
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Appendix A: Modelling the MIDI observations 

In this appendix we outline the method used to compare mod- 
els for geometry of the excess emission to the MIDI data. We 
consider symmetric Gaussians of varying FWHM and ring mod- 
els. The ring models are face-on axisymmetric discs with a fixed 
width and varying radius. The disc is assumed to be centred on 
the point-like star and has uniform surface brightness. The tested 
widths are dr/r = 0.2, 1 or 2, with different disc radii, parame- 
terised by r, the mid-point of the disc (so disc has inner radius 
r - dr/2, outer radius r + dr/2) with r between 5-200 mas. Due 
to the large range of possible non-axisymmetric disc models (for 
example discs with brightness asymmetries or clumps, or those 
inclined to the line of sight at all position angles) we consider 
only a restricted range of discs outside these face-on models. 
These were edge-on discs and those inclined at 45 and 70° to 
the line-of-sight. These discs were assumed to have a fixed ra- 
dius (80mas for HD69830, 90mas for 77 Corvi; sizes from SED 
fitting) and had position angles on sky of 0-180° East of North 
in steps of 15°. 

The excess region (Gaussian or ring) visibilities were deter- 
mined using the Interfere/metric Visibility Calculator available at 
http://www.mporzio.astro.it/%7Elicausi/IVC/. Determination of 
the total model visibility then follows from equation[3] assuming 
the visibility of the star is 1 at all wavelengths. We determined 
the visibility for the models across the wavelength ranges used 
in the characterisation of the dips (i.e. 8-9/im and 10-11.5/im, 
Tables |3] and |4]i. The models assume flux ratios for the star and 
disc components are accurately given by the IRS photometry and 
the Kurucz model profiles. 



A.1. Applying a goodness-of-fit test 

As the absolute levels of visibility are poorly constrained but the 
ratio of visibilities at different wavelengths is constrained much 
better (see sections 2.4, 3.1 and 3.2), we shall use comparison 
between the observed and model visibility ratios to determine 
the goodness-of-fit of different models. The sources of error for 



this comparison must be carefully considered to allow an accu- 
rate goodness-of-fit calculation. 

In this paper we have used the IRS photometry to calculate 
visibilities. Following equations[TJand[2] each observed visibility 
is calculated as 



V = /, 



corr,tar/ ^corr,cal 
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xK 
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When we consider the ratio of visibilities, we divide the two 
visibilities as averaged over the appropriate wavelength ranges 
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where the subscript 10 refers to the value of this term averaged 
over the wavelength range 10-1 1 .5/jm, and the subscript 8 refers 
to averaging over the range 8-9/jm. We now consider the errors 
from each term in the above equation. 

— , is from the correlated 



The first term in equation lA.2l j 
fluxes as measured in the fringe tracking. Errors in each wave- 
length range are given by variation in the five sub-integrations 
of the fringe tracking (as given in Tables [3] and |4]i and are prop- 
agated as o7, ratio // ratio = vV/.s/^s) 2 + (07,io//io) 2 - Th e second 
term, F ^ u " uo t is given by Kurucz model profile fits to the stan- 
dard stars. As the slope is considered to be well modelled by 
such profiles in the mid-infrared and there are no spectral fea- 
tures expected in this range for the standards we consider errors 
on this term to be negligible. The third term, /c " rr - c " u comes from 
the standard star correlated flux measurements. As shown in sec- 
tion 2.4 this ratio should always be 1 but has an error of 3.7%. In 
those instances in which an average calibration from two stan- 
dard star observations was used this error can be reduced by a 
factor of V2. The term p 01 '""' 8 is given by the IRS spectrum of 
each target. For this we must co nsider how well the IR S spectra 
captures the shape of a spectrum. lBeichman et alj d2006f) showed 
that for Sun-like stars with no excess emission deviations from 
a photospheric slope were 1.3% in the range 8.5-13/jm after re- 
moving the effects of absolute calibration error (by scaling to the 
photospheric model at 7-8//m). The stars included in this study 
are of the same or similar spectral type to HD69830 (which was 
in fact included in the sample), rj Corvi was just excluded by 
spectral type; the paper considered stars in the range F5-K5, but 
we do not consider an F2 star likely to be very different to this 
sample. We therefore use the error of 1 .3% as the error on the to- 
tal flux ratio for the science targets. The term is determined 

D Vcal,8 

by the assumed size of the standard star used for calibration and 
in these calculations is assumed to be perfect. In fact as can be 
seen in Table[TJthere is an uncertainty on the size of the standard 
stars. This uncertainty is very small however (typically O.Olmas) 
and potential error on the calibrator visibility ratio is thus very 
small (<0.01%). We therefore neglect this uncertainty in our er- 
ror calculations. 

To calculate the same ratio for each model we must use the 
following equation (c.f. equation[3): 



V, 



mod,ratio 



(-F*,tar,10/-Ftot,tar,l()) + (f disc,tar,lo/-Ftot,tar,l()) Vdisc,tar,10 
(*F*,tar,8/-F\ot,tar,8) + (^disc,tar,8/^tot,tar,8) Vdisc,tar,8 
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A.3) 



where it has been assumed the V* is 1 in both wavelength ranges 
(see section 3). The errors on each term for the model ratio 
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also require careful consideration. The first term in equation 
IA.3I , comes from Kurucz model profiles, and again as the 
stellar emission in this range can be considered to have a well- 
defined slope the error on this term is negligible. The second 
term, f ""•' ar - 8 , again comes from the IRS spectrum of this source, 

ftot,tar,10 

and thus we again consider the errors on this ratio to be 1.3% 
(see above). 

The errors on the final term in equation lA.3l are more com- 
plex. We consider that for each model the visibility Vdisc.tar is 
perfectly known. However errors on f tot, tar/ tar in both wave- 
length ranges come from a variety of sources, some of which are 
systematic across the different wavelength ranges considered. In 
order to determine the error on this factor we used a Monte- 
Carlo method. For each model (each with a fixed Vdisc,tar,8 and 
Vdi SC ,tar,io for each baseline) the value of this final factor was cal- 
culated over 1000 iterations, with the values of F lotMr and F + t ar 
chosen as follows. The values of F* it ar,8 and F*, t ar,io were deter- 
mined by using the Kurucz model profile for the correct spec- 
tral type (see Table[TJ). Scaling was determined by taking the^ 2 
minimisation over a range of scalings to find the best fit to the 
B- and V-band magnitudes of the stars as listed in the Hipparcos 
catalogue and J-, H- and K-band magnitudes from the 2MASS 
catalogue. The percentage points of the x 2 distribution were used 
to determine a l<x limit on this scaling. The scaling used in the 
Monte-Carlo testing was then taken from a normal distribution 
with a mean given by the best fitting scaling and standard de- 
viation at the level of the lcr error. The Fto^tar.s and F to t,tar,io 
values were determined by averaging over the appropriate wave- 
length ranges after the flux in each spectral channel was ran- 
domly sampled from a normal distribution with a mean of the 
IRS total flux in the channel and standard deviation the error in 
the channel. To approximate the effect of calibration error the 
values of F to t,tar,8 and f t ot,tar,io were then multiplied by a factor 
randomly sampled from a normal distribution with mean 1 and 
standard deviation 0.05 (c alibration uncertai nty for IRS typically 
5% at short wavelengths; iHines et al.l l2005): the same "calibra- 
tion value" was used for F to t,tar,8 and F to t,tar,io- The mean value 
of this final factor in equation IA.3I was checked against the its 
value assuming no errors on Ftot.tar and F+. t ar and found to be the 
same. The error term on this final factor was then taken from the 
standard deviation over the 1000 iterations, and was found to be 
typically at the 1 % level for both sources averaged over all disc 
model and baseline combinations. 

To test how well each model reproduces the observed vis- 
ibility ratios we use the y 1 goodness-of-fit test. Typically this 
takes the form E^ =1 (D - M) 2 /cr 2 where N is the number of data 
points, D is some observed data with associated error cr and M 
is a model with no error. As both V bs,ratio and V mo d,ratio have er- 
rors associated with them, we use the ratio of the observed to 
model visibility ratios and compare them to 1 (as if the model is 
a good fit to the data then V mo d,ratio ~ Kbs.ratio)- The errors from 
both visibility ratios must be included in the value of cr. Our \ 2 
calculation then becomes 

2 /Kibs,ratio(0/^mod,ratio(0 — 1 \ 2 < N 

X -g( " ) (A.4, 

where i labels the observations of each source, and cr, is the er- 
ror on Vobs.ratio/ V mo d,ratio- Consideration of equations lA.2l and lA31 
shows that the F to t,tar term cancels when taking the ratio of the 
observed and model visibility calculations, and thus the 1.3% er- 
ror term does not explicitly appear in our^ 2 calculations. Errors 
on the IRS photometry do appear implicitly in the Monte-Carlo 



calculation of the error on the final term in equation lA.3l through 
the calibration error and the statistical errors. The final value of 
cr comes from addition in quadrature of all the remaining error 
terms outlined in the above paragraphs. We use the percentage 
points of the^ 2 distribution to determine the values of r that pro- 
vide a reasonable fit to the observed data for both sources. The 
results of this modelling are discussed in sections 5.1 and 5.2 
and shown in Figures [8] and [10] 



